Density of states governs light scattering in photonic crystals 
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We describe a smooth transition from (fully ordered) photonic crystal to (fully disordered) pho- 
tonic glass that enables us to make an accurate measurement of the scattering mean free path in 
nanostructured media and, in turn, establishes the dominant role of the density of states. We have 
found one order of magnitude chromatic variation in the scattering mean free path in photonic 
crystals for just ~ 3% shift around the band-gap (~ 27 nm in wavelength). 
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Contemporary photonic science is capable of address- 
ing fundamental questions at the basis of light-matter 
interactions, of which the role of the photon density of 
states, Viuj), in nanostructured media is one of the most 
intriguing. 

Artificially engineered materials allow the control of 
light transport through interference in the internal nanos- 
tructure, rather than on the refraction in the body 
boundaries, engendering new materials properties. Pho- 
tonic crystals, in which the dielectric constant is period- 
ically modulated, manipulate electromagnetic states and 
the available phase space and control fundamental as- 
pects of light-matter interaction like light emission [l| 
and light transport much like semiconductors con- 
trol electrons. Redistribution and inhibition of the emis- 
sion from photonic crystals was proven [3], but uncon- 
ventional light transport in partially-disordered photonic 
crystals has only been hinted at by pioneering experi- 
ments 0, [5|. Other topologies, hke random media 0, 
correlated disordered Q or fractal [1], employ the ape- 
riodic subwavelength dielectric nanostructure to achieve 
similar light control for transport [9] and random lasing 
emission [10|. 

Light scattering by weak topological disorder in a pho- 
tonic crystal and the interplay between order and disor- 
der has yet to be fully understood and explored. As an 
important step, the relation between scattering extintion 
and 'D{uj) has just been theoretically derived As 
pointed out by John [2] a dramatic change in light diffu- 
sion can occur for frequencies in or around the band- 
gap and eventually Anderson localisation of light can 
be reached, the photonic conductor becoming an insu- 
lator 1^. In the quest for light localization, the first 
experiments focused on fully random media [l3| , only re- 
cently transverse localization has been reported in two- 
dimensional crystals with disorder 14l |. 

Even far from the localization regime, the scattering 
properties of Bloch modes, the periodic electromagnetic 
modes of a photonic crystal, are expected to be pro- 
foundly different from the diffusive modes encountered in 
conventional random media. Pioneering experiments on 
coherent backscattering [l^, and diffuse light transport 



[id . |l7| in photonic crystals searched for signatures of 
Bloch-mode mediated scattering but have merely shown 
standard light diffusion. Moreover, the experiments have 
been interpreted using a model that assumes no photonic 
modal dispersion but rather a modified reflectivity at the 
system boundaries [l8| . 

In this letter we study the scattering mean free path, 
£s, the fundamental building block for any wave transport 
model, for the special case of photonic crystals with a 
controlled amount of disorder. We report experimental 
evidence of strong chromatic dispersion of ig from band- 
edge to band-gap, and values of up to ^ 100 — 500 /im, 
i.e. ~ 300 times the lattice parameter (a), an order of 
magnitude higher than previously reported [HHIT^. 

Single scattering events in a system with modified light 
modes and density of states, as in a photonic crystal, are 
very different from those occurring in vacuum due to: a) 
an increase of light-matter interaction and thus of scat- 
tering by defects, when 'D{lu) is increased at the vicinity 
of band-edges and b) a suppression of the scattering chan- 
nels, i.e. a increase of is in the band-gap, where "D^uj) is 
strongly reduced. 

The scattering strength can be studied via simulta- 
neous reflection and transmission measurements, when 
absorption is negligible (the absorption length £a ~ 10 
m [19]) and for energies below the onset of diffraction 
(a/A 1.12 ^ where A is the hght wavelength). We 
assume that scattering losses follow Lambert-Beer's law, 
i.e. that after a thickness L, a ballistic beam attenuates 
as I{L) — Iq exp{~L/£s)- The intensity balance can then 
be expressed as 



r(L) + i?(L) =exp(-i/4). 



(1) 



where T{L) and R{L) are the ballistic transmission and 
reflection as a function of sample thickness in a given 
direction. 

The samples are Polymethil-metacrilate (PMMA, re- 
fractive index, n = 1.49) self-assembled face centred cu- 
bic (fee) photonic crystals [2l[ with controlled density of 
intentionally added vacancies (from 0% to 40%) (see 
Fig. [T^ and [T]d). These vacancies are obtained upon 
removal of given fractions of constituting spheres from 
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random lattice positions. 
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FIG. 1: Top panel: SEM images of photonic crystal with 
0% (a) and 40% (b) vacancy doping, (c) and (d) reflection 
and transmission from the corresponding samples of the above 
panel. 

The scattering mean free path in a random medium 
can be shown to be is = 1/ PsCr, where a is the scatter- 
ing cross-section and ps is the scatterers number density 
[^. is not only a measure of the quality of a pho- 
tonic structure but also the basic length-scale of a more 
complex picture of multiple scattering and light diffusion. 
While quantities like the transport mean free path or the 
diffusion constant are meaningful only in the context of 
the diffusion approximation, has a full meaning in any 
microscopic picture, regardless of the transport regime. 

In our system, the degree of extrinsic disorder (see 
scanning electron microscopy (SEM) images in Fig. [1^ 
and[T]D) can be very precisely and uniformly tuned, while 
keeping the sample thickness controlled. This allows us 
to develop a setup to measure Lambert-Beer's law for 
photonic crystals. We used a Fourier transform spec- 
trometer coupled to a microscope allowing to probe the 
scattering in the (111) direction while illuminating a 
constant-thickness area with a spot of ^ 80 /im. The 
spectra (see figures [TJ; and[TJl) are taken in adjacent re- 
gions which are visible by optical microscope inspection 
as terraces on the sample surface. The thickness of such 
layers is precisely measured with an uncertainty of 2%, 
via the density of Fabry-Perot fringes that occur for the 
interference of the light reflected from the front and rear 
faces of the sample. 

Figure [2^ shows the measured ln(r + R) for three dif- 
ferent degrees of vacancy doping Q i.e. for different de- 
grees of extrinsic disorder, at a wavelength of 633 nm and 
for spheres of 237 nm in diameter (a/ A — 0.52). In this 
type of representation, the slope yields directly (— £s)~^ 
according to Eq. [T] This wavelength is chosen to exem- 



plify a spectral region where no photonic band features 
are present, as, at such a low energy, the photonic crystal 
band dispersion is the same as in a uniform homogeneous 
effective medium. In figure [5^ three scattering regimes 
are clearly distinguishable. For thicknesses lower than 
ca. 10 layers, (regime I), up to 25 — 30% of the inci- 
dent light is scattered due to surface effects, in the form 
of stacking patterns [l^l , high lattice displacements and, 
even stacking order arrangements [1^. When the sec- 
ond regime (II) sets in, the slope reaches a station- 
ary value that characterises the photonic crystal. Eq. [1] 
holds and scattering losses scale with sample thickness 
like ~ ey.'p{L/tg). Finally, for larger thicknesses, a third 
scattering regime (III) appears in Fig. [5^. Apparently, 
for thick samples > 50 layers, the self-assembling process 
loses its effectiveness, as it evident by the increase of in- 
trinsic disorder and by the cracks that appear on thick 
sample as they start to lift from the substrate. 

The physical picture we propose can be checked against 
consistency if additional disorder is added to the pho- 
tonic crystals. This can be done by doping the original 
photonic crystal with a controlled concentration of va- 
cancies. At a wavelength of 633 nm the values of Ig, cal- 
culated from the fit to the Lambert-Beer law, are plotted 
as a function of vacancy density in Fig. ^Bp. The "per- 
fect" crystal (that with 0% added vacancies) is highly 
ordered as it presents a scattering mean free path of 63 
/im, hundreds of times the lattice constant (in this case 
a ~ 0.33 /im), and in particular much larger than the 
Bragg length (of the thick sample) (l^l that in our case 
is Lb — (3.8 ± 0.3) /im. An addition of a very little 
amount of defects rapidly decreases the mean free path. 
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FIG. 2; (a) Plot of ln{R + T) as a, function of the sample 
thickness, at A = 633 nm, for different vacancy density doped 
photonic crystals (from 0% to 40% vacancies doping), of 237 
nm diameter, (b) is obtained from linear fitting of the slope, 
for regime II. It also shows the Bragg length {Lb) in the case 
of pv — 0% as shaded area. 
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hence the quaUty of the crystal. In this figure, the inverse 
scattering mean free path scales linearly with the vacancy 
concentration py, as shown by the black line, which is a 
fit for ij^ = Po<tq + pyOy where po and p^ are the density 
of intrinsic and intentionally added scatterers and (Tq a-nd 
Cy their scattering cross-section respectively. From the 
fit of £s{pv) as a function of the vacancy concentration 
we can estimate ay = (0.016 ± 0.002) fim^. 
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FIG. 3: Figure shows is as a function of the light wavelength 
for 0%, 10% and 40% vacancy doped photonic crystals with 
d = 237 nm. The position of the pseudo-gap is shaded in 
cyan. The dotted line shows the dependence of^s far 
from the band-gap. 

Fig. [3] shows the strong chromatic dispersion of £s{^) 
in the visible range. This is the signature of the photonic 
crystal. In the low energy side of the pseudo-bandgap, 
£s takes on a value of the order of ~ 100 pm for sphere 
diameter d = 237 nm and ~ 500 pm for d ~ 600 nm (not 
shown here), the largest values reported so far. Previ- 
ous experiments 13, 3 12 1 have measured it, the trans- 
port mean free path, from very thick ('^ 200 pm) pho- 
tonic crystals grown by natural sedimentation 3 12 1 or 
centrifugation [15], and found values in the range 7 pui 
< it < 20/j,m. is is in general smaller than it and, there- 
fore, the values found in our experiment represent a much 
higher degree of ordering than previous reports. Far from 
the band-gap, is{i^) varies as ~ dependence that 

has been confirmed also in previous experiments [l5| and 
attributed to Rayleig-Gans type of scattering. 

A simple model for a point-like scatterer in a photonic 
crystal can explain the chromatic dispersion of is- We as- 
sume that the scattering is isotropic and, therefore, that 
the imaginary part of the green tensor G(r, r') can be 
approximated by Im{G(r,r')} ~ {ttc^ /2uj)'D{r,uj)I. We 
denote the k contribution to the local density of states (or 
"projected" density of states) by X'k(r, w), and calculat- 
ing the power radiated by a dipole we obtain a (spatially 
averaged) scattering cross section 

ak(cc>) ~ F{lj)uj^V{u;)Vi,{lj), (2) 

where F^lo) is a form factor which takes into account cor- 
rections beyond Rayleigh scattering and k the wavevec- 



tor in the incident direction. In our simple model, the 
form factor can be replaced by a Rayleigh-Gans factor 
F{llj) Ld^"^ and the polarizability of the scatterer can 
be considered independent of frequency [13. Eq. [2]states 
that the scattering cross-section, and hence the scatter- 
ing mean free path, has a strong dependence on the total 
2?(a;) and projected 'D\^{lo). This dependence typically 
disappears in ordinary random media for which the pho- 
tonic modes are isotropic and energetically smooth, but is 
very important for photonic crystals. At the band-edges 
of our photonic crystals, is{^) has a sharp decrease of a 
factor of up 4 to is{a/\ = 0.59) = 11 ± 1 pm and then 
it shoots up almost an order of magnitude in the band- 
gap to is{a/X = 0.62) = 81 ± 40 pu\. Such an 8-fold 
increase occurs within just 0.03 in a/A and ^^27 nm in 
wavelength, around the photonic band-gap. Again, as 
a comparison, we show in Fig. [3] the frequency depen- 
dence of is for the 10% and the 40% vacancy photonic 
crystal, the latter can be considered fully disordered. As 
the vacancy doping is increased, the profile is smoothed. 
Firstly the band-edge effect on is disappears, as these 
standing-wave-like states are very sensitive to disorder. 
Then, for the 40% vacancy case, the effects on V{Ld) are 
washed out and the only feature in is occurs at the posi- 
tion of the first Mie resonance of the individual dielectric 
spheres ^19;] . This weak energy dependence of is is likely 
to be the only residual effect in a very disordered opal, as 
tho se g rown by natural sedimentation or centrifugation 
(Til Ho . \y\ , which present superficial iridescence but are 
largely bulk-disordered and exhibit standard light diffu- 
sion. 
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FIG. 4: (a) Inverse of the total "Diuj) and the projected 
^'fc=(i,i,i) ('^) a-long the incident direction (F — L direction) 
are plotted, (b) Enhancement factor, £s{a/X) /ls{a/\ = 
0.4), for two opals with no vacancy doping and the prod- 
uct CDkiuj) — 1). The position of the pseudo-gap is shaded in 
cyan. 



Figure 2^ shows the inverse of the total density of 
states (dark blue dotted line) that has a very weak mod- 
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ulation at the gap together with ^^^{uj) for propagation 
parallel to k = (111) that does have a strong variation 
at the gap (violet full line). Pk(^) is calculated from the 
inverse of the group velocity [25| ?;g(111) from ref. [1^. 

Figure IHd shows the enhancement factor defined as the 
ratio of £s to its value far from the gap, €s(a/A)/£s(a/A — 
0.4), for two different samples composed by PMMA 
spheres of 237 nm (black curve) and 600 nm (green curve) 
in diameter respectively. The enhancement factor points 
out the existence of a photonic pseudo-gap and reveals 
the variation of the density of states in the photonic crys- 
tal. A clear resonant behavior is evident. The variation 
in is is 8-fold for d = 237 nm and 20-fold for d = 600 
nm, which we attribute to the superior quality of the 
lattice. In figure 4b we plot also 1/P(w)2?(u;)k (dashed 
pink curve) that, from Eq. [21 is expected to reproduce 
the shape of energy dependence of £s{^^)- A fair agree- 
ment between theory and experiment is obtained, and 
the qualitative behavior is well captured by our simple 
model. Although both I7(a;) and X'k(w) contribute to 
the strong variation of is , it is evident that the principal 
factor responsible for a change in is (111) is X'k('jj)- 

Our simple and qualitative model accounts remarkably 
well for the shape of the measured is{uj) although it does 
not account for the asymmetry of its dispersion in the 
photonic gap. This effect is related to the available scat- 
tering states in other crystallographic directions close to 
the incident one 27|. Our measurement indicates the 



need for a more complete theoretical model that should 
account for all the modes in all directions with the right 
scattering probability. 

An increase of is in the band-gap and a decrease at 
the band-edge reflects the modified phase space avail- 
able Ak for light scattering when the photonic modes 
are concentrated around few fc-directions or the available 
scattered states reduced. This is consistent with John's 
seminal prediction of a need for a modified loffe-Regel cri- 
terion fe?! for scattering in photonic crystals, to include 
Afc. In addition, here we show that as the phase-space is 
modified, also is is altered: light scattering in photonic 
crystals is richer than in conventional amorphous media. 
Complete photonic band-gap materials, like Si inverted 
opals, would amplify the effect here presented and could 
be proper candidates to observe Anderson localisation of 
light. 

In conclusion, we show that a controlled smooth tran- 
sition from ballistic to diffuse transport in photonic crys- 
tals can be induced by the introduction of extrinsic dis- 
order. We find that the strength of scattering is closely 
related to the density of states, which induces immense, 
up to 20-fold, variations in the scattering mean free path. 
We propose is as a robust, easy to measure, figure of 
merit in assessing the quality of photonic crystals for 
technological applications. The possibility of controlling 
light scattering and diffusion in nano-structured optical 



media has important implications not only to test the 
quality of photonic devices, but also to properly address 
the proximity to the onset of Anderson localization in 
disordered lattices, or for the spectral control of lasing 
emission from disordered/ordered active media [lo| . 
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